Magnetoconductance through a vibrating molecule in the Kondo regime 
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The effect of a magnetic field on the equilibrium spectral and transport properties of a single- 
molecule junction is studied using the numerical renormalization group method. The molecule is 
described by the Anderson-Holstein model in which a single vibrational mode is coupled to the 
electron density. The effect of an applied magnetic field on the conductance in the Kondo regime 
is qualitatively different in the weak and strong electron-phonon coupling regimes. In the former 
case, the Kondo resonance is split and the conductance is strongly suppressed by a magnetic field 
gUsB > ksTK, with Tk the Kondo temperature. In the strong electron-phonon coupling regime a 
charge analog of the Kondo effect develops. In this case the Kondo resonance is not split by the 
field and the conductance in the Kondo regime is enhanced in a broad range of values of B. 
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I. INTRODUCTION 

Transport through molecular nanodevicesi has re- 
cently been a subject of intense experimental^^iSiS*^ 
and theoretical ll ' ) i 11 ' 12 i 1 ?' 14 i 1 5' 1 9i 1 7 investigation. These 
systems are not only of great interest because of their po- 
tential for technological applications but they also pose 
some interesting an challenging theoretical problems. 
They arise because in molecular nanodevices electron- 
electron (e-e) and electron-phonon {e-ph) interactions 
are present simultaneously and their interplay may lead 
to new physics if their strengths are comparable. 

Effects of e-e interactions in single-molecule devices 
were seen through the observation of Coulomb blockade 
peaks in the conductance. At these peaks the charge of 
the molecule changes by one electron while it stays well 
defined in the valleys between them£ 

The Kondo effect was also observed in Ceo sin- 
gle molecules coupled to metallic^ and ferromagnetic 
electrodes^ It occurs in the valleys in which the charge 
of the molecule is close to an odd integer. Its experimen- 
tal manifestation is an increase of the conductance with 
decreasing temperature and a zero-bias peak in the dif- 
ferential conductance below a characteristic temperature 
scale, the Kondo temperature Tk. 

Furthermore, some features in the I-V characteristics 
of a vibrating Cqq single-molecule transistor were inter- 
preted as being a manifestation of the effects of the e-ph 
interaction^ 

It was shown in previous worki^i^ that the inter- 
play of e-e and e-ph interactions has important effects 
on the physical properties of the system in the Kondo 
regime. There are two different cases characterized by 
the strength of the e-ph coupling. When the latter is 
weak (in a sense that will be defined below), the ground 
state of the isolated molecule in the odd electron num- 
ber valleys is a spin doublet. Then, the conventional 
Kondo effect appears when the molecule is coupled to 
metallic leads. The presence of the e-ph interaction re- 
sults in an increase of the spin fluctuation energy and of 
the amplitude of charge fluctuations. In the strong e-ph 



coupling regime, polaronic effects result in an effective 
attractive interaction U e s between the electrons and this 
leads to qualitatively new physics. If the charge in the 
device is close to an odd integer the ground state of the 
isolated molecule is a spinless charge doublet. This is 
the regime in which the charge Kondo effect develops. 
It differs from the conventional Kondo effect in that the 
low lying excitations of the system are charge fluctua- 
tions whereas spin fluctuations become increasingly stiff 
as the e-ph coupling increases. The charge Kondo effect 
enhances the transport through the molecule and leads 
to the appearance of a narrow peak in the conductance 
function of the gate voltage. 

In this paper we study the properties of the device in 
an applied magnetic field. We found that these are quite 
different at weak and strong e-ph coupling. At weak cou- 
pling the system behaves qualitatively as in the absence 
of the latter iiS The magnetic field breaks the spin symme- 
try of the molecule's ground state and splits the Kondo 
peak in the density of states thus suppressing the conduc- 
tance for fields g^isB > fcsTK. At strong coupling, the 
molecule's low lying states are spinless and are thus insen- 
sitive to the magnetic field. The Kondo peak is not split 
by the latter and there is almost perfect transmission at 
low temperatures for fields B < B cl i t = \U c s \ /g^B- Fur- 
thermore, within this range of field values, the width of 
the conductance vs. gate voltage peak increases rapidly 
with B. Beyond B cvlt there is level crossing and the con- 
ductance peak splits into two peaks that are strongly spin 
polarized. 

The rest of this paper is organized as follows. In Sec- 
tion [n] we present the model Hamiltonian and some gen- 
eral considerations that follow from exact Fermi liquid 
relationships at T = 0. In Section ITTT1 we present Numer- 
ical Renormalization Group results for the spin resolved 
spectral densities and conductance in the presence of a lo- 
cal magnetic field at finite temperature and as a function 
of gate voltage. A summary and our concluding remarks 
are presented in Section ITVl 
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II. MODEL HAMILTONIAN 

We consider a model of a molecule with a single rel- 
evant electronic level coupled to left (L) and right (R) 
metallic electrodes and a vibrational mode of frequency 
ojo linearly coupled to the charge fluctuations of the 
molecule. The Hamiltonian of the system is 

H = Hm + He + Hm-e + Hmf, (1) 

where the first three terms describe the isolated molecule, 
the electrodes and their coupling, respectively. The last 
term is the coupling of the spin of the molecule to an 
external magnetic field acting only on the molecule. We 
have 

H M = s d hd + Un d ^h dl - A (n d - 1) (a + a f ) 

+ ujna^a, (2) 
He = £ <?»( k ) c La c k CTQ : (3) 

k,cr,Q={L,R} 

,(4) 

k,<r,a={L,R} 

Hmf = gfJ-BBS^. (5) 

Here, fi d — ^<rd\d a and S% — |(d|rff — d\d-i) are 
the charge and spin operators of the molecule, respec- 
tively. The position of the electronic molecular level rel- 
ative to the Fermi level of the electrodes is denoted by 
e d and U is the Coulomb repulsion between two elec- 
trons that occupy the same molecular level. Finally, cr 
is the creation operator of a phonon of frequency wo and 
A is the e-ph coupling constant. We consider for sim- 
plicity the case of identical electrodes with dispersion 
er(1c) = Eij(k) = e(k) and a k-independent tunnelling 
matrix element. We further assume that the conduction 
band is symmetric around the chemical potential and de- 
note its half- width by D. In the following we use units 
such that fcs = D = gfj,B = 1- 

The eigenstates of the Hamiltonian Hm + Hmf of the 
isolated molecule are direct products of electronic states 
(denoted by a subscript e) and oscillator states. The 
eigenf unctions and their energies can be written explic- 
itly: 

|0,m) = tf-|0) o |m>, E^-g + ma/o, 

|cr, m) = |cr) e \m), E m = e d + oB + mwo, (6) 

|2, m) = U+ |U) e \m), El = -±i + 2e d + U + mw , 

where U — exp [±A/cjo (a — a)] , \m) is the m-th ex- 
cited state of the harmonic oscillator and a — ±1/2. The 
ground state energies of the empty and doubly occupied 
electronic states are reduced by the polaronic shift A 2 /loq. 

At zero field the differences between the occupied 
and empty states E m — E m — E m and E m — E m are 
E? n = £d + A 2 /cl> and Ef n — 2e d + U, respectively. The 
effective interaction between two electrons occupying the 
same state is thus E m — 2E m = U — 2A 2 /w = U e g. 



The sign of this quantity determines the low energy 
physics of the system and its response to an applied 
field. For positive U c g, the ground state of the isolated 
molecule is a spin doublet and the Kondo effect occurs 
at low temperatures when the molecule is coupled to the 
leads. For negative U e g the ground state of the isolated 
molecule is either empty, doubly occupied or a charge 
doublet. In the latter case the coupling of the molecule 
to the leads generates the charge Kondo effect with an 
anisotropic Kondo (AK) exchange couplingiiSiiSiSli 

In this paper we only study the case of paramagnetic 
contacts and, therefore, we set F ct l(w) = T^^u)) = 
|r(w), where T(u) = 2n £ fe 5 {oj - e(k)) Vg. The spin 
resolved conductance of the molecular junction at zero 
bias ia 21 i 22 i 2 'V 4 i 25 

G a = j ir J™ dco f-^^) fe(w)rM, (7) 

where p d(J {u) = — ir~ 1 ImQ da -(u>), Q da {uj) is the retarded 
electronic Green function of the molecule in the presence 
of the leads and f{ui) is the Fermi distribution. 

At zero temperature, in the wideband limit, the di- 
mensionless conductance satisfies the Fermi liquid rela- 
tionship 

9 = = ^y^T £W(0) = V sin 2 (nn da ) , (8) 

(jro 

cr cr 

where Go = e 2 /h is the quantum of conductance, n d(T = 
(h d<7 ), T = r(0), and the second equality follows from 
Luttinger's theorem^ generalized to the case in which 
e-ph interactions are present^ 

Consider first the case ^ CT n da — 1 corresponding to 
the electron- hole symmetric situation e d — —U/2. The 
conductance in zero field is then g = 2. In the presence of 
the applied field B the molecule acquires a magnetization 
m d = \ (n>dl ~ n d|) and Eq. © becomes 

g = ^sin 2 [7r(l/2 - 2am d )\ = 2 cos 2 (irm d ) . (9) 

In the low field limit we have, 

5 «2(l-7T 2 m 2 )«2(l-^ 2 X ^ 2 ), (10) 

were Xd = dm d /dB is the local spin susceptibility. In 
the large field limit we have m d — > 1/2 and the conduc- 
tance vanishes. The characteristic field beyond which the 
conductance is suppressed is B* ~ X d ■ 

For U c s large and positive x d X ~ Tk and increases 
slowly with increasing A>i£*i& Therefore, B* ~ Tk is small 
in the limit U g 3> T. For B > B* the conductance at 
£d ~ — U/2 is thus suppressed and g acquires a double- 
peak structure as a function of the gate voltage. The zero 
temperature conductance at the positions of the peaks is 

For U e g large and negative x d 1 ~ I U e & \ and B* is thus 
largeii^i^ In this case, a single conductance peak with 
g = 2 is expected in a broad range of magnetic fields. 
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FIG. 1: Total and spin-resolved spectral densities at T = 
for U e ff > 0. Parameters values are Ed = —0.04, Y — 0.00393, 
U = 0.1, U) = 0.02, and U cB = 0.0375. T K = 8 x 10" 5 . Full 
line: total spectral density for B = 0. Dotted and dashed 
lines: spin down and spin up components, respectively, for 
B ~ Tk. Inset: structure of the Kondo peak. Full thick 
and thin lines: total spectral density for B — and B ~ Tk, 
respectively. Dotted and dashed lines: majority and minority 
spin components of the spectral density, respectively. 

In zero field but off electron-hole symmetry, Ed = 
— U /2 + 8sd, we have: 

J»2(l- 7T 2 Sn 2 d /4) w 2 (1 - 7T 2 x^/4), (11) 




FIG. 2: Zero temperature spectral density in the region of 
the Kondo peak and conductance as a function of tempera- 
ture and magnetic field at the electron- hole symmetric point. 
In all cases Ed = —0.05, V = 0.01, and ujq = 0.05. (a): 
Zero temperature spectral density for U e g — 0.051 and three 
values of the magnetic field. Tk = 0.004. (b) Zero tem- 
perature spectral density for U c s = —0.044 and three val- 
ues of the magnetic field. Tak = 1-1 X 10~ 4 . (c) Conduc- 
tance as a function of the magnetic field and temperature 
for UeB — 0.051. Temperatures are T — (top curve) and 
T = 0.006 x 2~ N ,N = 0,1,2,..., 5 (from bottom to top), 
(c) Conductance as a function of the magnetic field and tem- 
perature for U e B — —0.044. Temperatures are as in (b) with 
AT = 0,1,2,...,9. 



where Xc = —dn/d£d is the charge susceptibility and 
$ n d = ^2rj n da — 1- The width 5e* d of the conductance 
peak is ~ Xc 1 ~ ^ak, the Kondo temperature of the ef- 
fective anisotropic Kondo model that is small for very 



large |J7 e ff|, 5e* d ~ exp 



M In the 



(ir\U cS \* /ATujo) 
same limit the field-dependence of 5e* d can also be eval- 
uated analytically (see Eqs. i|13|) an d l|14|) below). We 
find 



KM 



exp [ttB 2 / (16rw )] 



(12) 



The width of the conductance peak thus increases expo- 
nentially with B 2 in the strong coupling limit. 



III. NUMERICAL RESULTS 

In this section we present numerical results for the spin 
resolved spectral density Pda{u) = — 7r _1 Im Qd a (u>) and 
conductance obtained using the Numerical Rcnormaliza- 
tion Group method (NRG)2£i2I modified to include e-ph 
couplingSS and to calculate accuratel y th e spectral den- 
sity Pd<j(w) using the method of Ref. 29]. 

Figure ^ represents the total spectral density in zero 
field (full line) and its spin resolved components (dotted 
and dashed lines) in the presence of a small magnetic 



field B ~ Tk for a case with U g > 0. Parameters are 
e d = -0.04, T = 0.00393, U = 0.1, U eS = 0.0375 and 
cjo = 0.02. The system is off the electron-hole symmetric 
point and Sed > 0, rid < 1< The density of states is thus 
asymmetric about the Fermi level, Pd(w) ^ Pd{— w). The 
total spectral function is represented by the full line. We 
observe the presence a very narrow Kondo peak near to = 
and additional structure on either side of this feature. 
The Kondo temperature, determined from the full-width 
at half maximum of the central feature, is Tk ~8 x 10~ D 
for this set of parameters. The structure on either side of 
the Kondo peak are phonon side bands as can be seen by 
analyzing the pole structure of G^, the retarded Green 
function in the atomic limit (i.e., at r = 0)i2& 

The vertical lines in Fig. ^ represent the positions of 
the poles of G%j given in Eq. (15) of Reference 13 It 
can be seen that they are in fairly good agreement with 
the positions of the peaks observed away from the Kondo 
feature in Fig. ^ 

The presence of a field B^Tk has drastic effects on the 
spin resolved spectral spectral function. Due to the high 
magnetic susceptibility of the molecule (oc 1/Tk) the lat- 
ter acquires a large magnetization. Most of the electron 
density is composed of majority spins (down spins) and 
the population of the minority spins (up spins) is small. 
The spectral function of the majority spins (represented 
by dotted lines) is thus large for oj < whereas that of 
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the minority spins (represented by dashed lines) is large 
for u>0. 

The inset to the figure shows in detail the effect of the 
magnetic field on the Kondo peak. The thick and thin full 
lines represent the total spectral function in the central 
region in zero field and for _B^Tk, respectively. It can 
be seen that the Kondo peak splits and its amplitude is 
reduced anticipating the suppression of the Kondo effect 
at higher fields. The dotted and dashed lines represent 
the spin down and spin up components of pd(u>), respec- 
tively. It is seen that position of the two Kondo peaks 
shift in opposite directions by an amount ~ ±B. 19 

Figure [3 displays the spin resolved spectral density at 
T = for three values of the magnetic field (a,b), and 
the conductance as a function of magnetic field and tem- 
perature (c,d) at the electron-hole symmetric point. The 
parameters are T = 0.01, U — 0.1, and ujq = 0.05. We 
only show the minority-spin spectral density pd\ since 
electron-hole symmetry and spin reversal invariance lead 
to the relationship pdi{w) = Pd]{—u). 

It can be seen that the Kondo peak behaves differ- 
ently in a magnetic field in the cases shown in Fig. Ufa) 
(U cS = 0.051 > 0) and Fig. Hb) (U eS = -0.044 < 0). 
Whereas in the former case the Kondo peak shifts and 
its hight decreases upon applying a small magnetic field, 
in the latter case, for fields B < \U e s\, its position and 
height remain unchanged. Note that the shape of the 
peak is asymmetric and its width increases sharply upon 
applying the field. The spectral asymmetry is due to the 
decrease of the density of minority spins in the presence 
of the field as discussed above. The field-dependence of 
the width of the Kondo peak will be discussed below. 

For even larger fields, B> \U e s\, there is a level cross- 
ing and the charge doublet is no longer the ground state 
of the isolated molecule. The charge Kondo effect can 
not take place and the associated peak in the density of 
states disappears. 

The dependence of the conductance on field and tem- 
perature for the cases of positive and negative U c s is 
shown in Figs.|2c) and^d), respectively. 

For positive U e g, the field- induced suppression of spec- 
tral weight at the Fermi level seen in Fig. |2Ja) translates 
into a suppression of the zero temperature conductance 
as can be seen in Figl^c). As a function of field this 
quantity exhibits a peak of width ~ Xk 7 as expected 
from Eq. i|10[l • With increasing temperature the conduc- 
tance decreases further. 

The case of negative U e g is shown in Fig. Hfd). At tem- 
peratures much lower than the zero-field value of Tak we 
observe a wide plateau of width ~ 2|J7 e ff|- The region of 
the plateau is that in which charge Kondo effect exists 
as discussed above and this leads to a high conductance. 
With increasing temperature the conductance first de- 
creases in the central region where Tak is the lowest (see 
below). 

At high temperatures, two peaks are seen at B ~ 
±|[/ e ff |. These peaks are associated to the crossing of the 
levels of the isolated molecule [cf. Eq. ijfJJl]: for \B \ > \ U c a\ 




SE d /T AK & d /T AK 

FIG. 3: Zero temperature conductance G versus gate voltage 
Ed for different values of the magnetic field B in the strong 
e-ph coupling regime (a-d). Solid lines: total conductance. 
Dashed and dotted lines: majority (down) and minority (up) 
spins, respectively. Other parameters as in Fig|5Jb). 

the ground state of the molecule is either ||,0) or |T,0), 
depending on the sign of the field. For \B\ < \U c g\, in- 
stead, the ground state of the molecule is spinless. These 
peaks are thus the analogs of the familiar Coulomb block- 
ade peaks that appear for positive U e s upon varying the 
charge. For negative U e g and in the electron-hole sym- 
metric situation that we are discussing the charge is fixed, 
however, and it is the spin of the molecule that changes at 
the peaks by an amount AS Z — 1/2 when the magnetic 
field varies. 

We now discuss the spin resolved conductance at T = 
as a function of gate voltage (i.e., of Ed) for the case 
Ucs < and several values of the magnetic field. Fig- 
ure E[a) shows the conductance at zero field. When the 
molecular level departs from the symmetric position, the 
molecule charge-polarizes and the conductance decreases 
according to Eq. <| 1 1 f) . As a consequence, the conduc- 
tance exhibits a narrow peak of width ~ Tak- 

When a small magnetic field is applied, Tak increases 
as mentioned above and the peak in the conductance 
broadens accordingly [cf. Fig. EJb)]. There is also a 
spin polarization of the current that changes sign at the 
symmetric point. This can be understood as follows: for 
negative but small Sed the charge of the molecule in- 
creases and the current is dominated by charge fluctua- 
tions between the states |2,0) and ||,0). This leads to 
a predominantly spin-up current. For positive but small 
Ssd the charge in the molecule decreases and the relevant 
fluctuations are now those between |0, 0) and ||,0) and 
the spin-polarization of the current is reversed. 

The width of the conductance peak further increases 
with the magnetic field but the peak height remains re- 
mains nearly unchanged as seen in Fig. Etc) . Notice the 
strong spin-polarization of the current reflected in the 
strong asymmetry of the peaks in the spin resolved con- 
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FIG. 4: Magnetic field dependence of Tak computed from 
Eqs. (2) and for the same values of the parameters as 
in Fig. [3] The thin line is a linear fit. 



ductance. In this regime the device acts as a very effective 



spin filter and the spin polarization can be reversed by 
a small change in the gate voltage. For magnetic fields 
larger than the spin gap B > \U g\, the ground state of 
the isolated molecule changes as mentioned above and 
the total conductance shows a two-peak structure. The 
heights of these are reduced by half as shown in Fig.EJd). 
These peaks are located at dej, ~ ± (B — \U c g\) /2 where 
there is a level crossing of the molecular states. The 
conductance peaks are almost fully spin-polarized in this 
limit. 

We conclude by discussing in detail the field- 
dependence of the charge Kondo temperature Tak- As 
shown in Reference ll6l the low-energy effective model for 
the charge fluctuations in the strong e-ph phonon cou- 
pling case is the anisotropic Kondo model with couplings 
Jii and J± and Kondo temperature 



,30,31 



Tak ~ (Jx/J\\)^ 



(13) 



In the present case, in the presence of a magnetic field, 
the coupling constants arei& 



J 



Jo 
2 



^ 2\ 2 /u U 

.,o=±l/2 1 1 



(A/o; ) 



(A/w ) 2m /ml 



1 -2aB/U + 2muj /U 



J± 



Jo 
2 



(-1) 



^ 2\ 2 /uj a U 

,,o=±l/2 ' 



(X/ujq) 2 " 1 /ml 



l-2aB/U + 2mw /[/' 



(14) 



r 



where Jo = 8 T/(irUpo). 

For very large e-ph coupling the expressions above can 
be evaluated analytically leading to the result quoted in 
Eq. (|12fl in Sec. [H] In the numerical calculations reported 
here we used moderate values of A and the expressions for 
the couplings in Eq. (|14|) must be evaluated numerically. 

Figure ^represents the field dependence of Tak for the 
same parameters as in Fig. [21 It can be seen that in this 
regime Tak still increases exponentially with B. 



IV. CONCLUSIONS AND SUMMARY 

We analyzed the equilibrium transport and spectral 
properties of a model molecular transistor when the spin 
symmetry is broken by applying an external magnetic 
field. For weak e-ph coupling, in zero field, the conduc- 
tance shows the well known Coulomb blockade peaks at 
high temperature and the enhancement of the conduc- 
tance due to the Kondo effect at low temperatures. An 
external magnetic field B > Tk suppresses the conduc- 
tance as a result of the splitting of the Kondo peak. The 
properties of the conductance at zero bias are qualita- 
tively the same as those found in the case in which there 
is no e-ph coupling but with a set of phonon renormal- 
ized parameters. The spectral density, however, presents 



a set of additional peaks separated by the phonon energy 
ujq, that can be understood analyzing the level struc- 
ture of the isolated molecule. The properties of the spec- 
tral density are experimentally accessible in transport ex- 
periments at finite but small source-drain voltages or in 
scanning tunneling microscopy (STM) experiments for 
molecules on surfaces. 

In the strong e-ph coupling regime the behavior of 
the system is qualitatively different. The odd Coulomb 
blockade valleys are completely suppressed and the low 
temperature physics is dominated by charge fluctuations 
instead of spin fluctuations. In a small range of gate volt- 
ages the charge Kondo effect takes place and results in 
an enhancement of the conductance at low temperatures. 
In contrast with the previous case, the Kondo peak is 
not split by a magnetic field. Instead, for a broad range 
of magnetic fields, it broadens as the latter increases. 
For large enough fields or temperatures a spin blockade 
regime is observed. 

The magnetic field dependence of the conductance of 
a molecular transistor with strong e-ph coupling resem- 
bles in many ways the gate voltage dependence of the 
conductance of a device with weak e-ph coupling in zero 
field. The two cases are related to each other through the 
transformation (S5) — > (fid) and B — > V g . This transfor- 
mation is not an exact one because the symmetries in the 
spin and charge sectors of the Hilbert space are different 
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but it gives a good qualitative physical picture. 
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